Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions.

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output.

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Perturbations in 3D genome organization can promote acquired drug resistance.

Acquired drug resistance is a major problem in the treatment of cancer. hTERT-immortalized, untransformed RPE-1 cells can acquire resistance to Taxol by derepressing the ABCB1 gene, encoding for the multidrug transporter P-gP. Here, we investigate how the ABCB1 gene is derepressed. ABCB1 activation is associated with reduced H3K9 trimethylation, increased H3K27 acetylation, and ABCB1 displacement from the nuclear lamina. While altering DNA methylation and H3K27 methylation had no major impact on ABCB1 expression, nor did it promote resistance, disrupting the nuclear lamina component Lamin B Receptor did promote the acquisition of a Taxol-resistant phenotype in a subset of cells. CRISPRa-mediated gene activation supported the notion that lamina dissociation influences ABCB1 derepression. We propose a model in which nuclear lamina dissociation of a repressed gene allows for its activation, implying that deregulation of the 3D genome topology could play an important role in tumor evolution and the acquisition of drug resistance.

Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential.

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients' outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.

MYC promotes immune-suppression in triple-negative breast cancer via inhibition of interferon signaling.

The limited efficacy of immune checkpoint inhibitor treatment in triple-negative breast cancer (TNBC) patients is attributed to sparse or unresponsive tumor-infiltrating lymphocytes, but the mechanisms that lead to a therapy resistant tumor immune microenvironment are incompletely known. Here we show a strong correlation between MYC expression and loss of immune signatures in human TNBC. In mouse models of TNBC proficient or deficient of breast cancer type 1 susceptibility gene (BRCA1), MYC overexpression dramatically decreases lymphocyte infiltration in tumors, along with immune signature remodelling. MYC-mediated suppression of inflammatory signalling induced by BRCA1/2 inactivation is confirmed in human TNBC cell lines. Moreover, MYC overexpression prevents the recruitment and activation of lymphocytes in both human and mouse TNBC co-culture models. Chromatin-immunoprecipitation-sequencing reveals that MYC, together with its co-repressor MIZ1, directly binds promoters of multiple interferon-signalling genes, resulting in their downregulation. MYC overexpression thus counters tumor growth inhibition by a Stimulator of Interferon Genes (STING) agonist via suppressing induction of interferon signalling. Together, our data reveal that MYC suppresses innate immunity and facilitates tumor immune escape, explaining the poor immunogenicity of MYC-overexpressing TNBCs.

Double-strand break toxicity is chromatin context independent.

Cells respond to double-strand breaks (DSBs) by activating DNA damage response pathways, including cell cycle arrest. We have previously shown that a single double-strand break generated via CRISPR/Cas9 is sufficient to delay cell cycle progression and compromise cell viability. However, we also found that the cellular response to DSBs can vary, independent of the number of lesions. This implies that not all DSBs are equally toxic, and raises the question if the location of a single double-strand break could influence its toxicity. To systematically investigate if DSB-location is a determinant of toxicity we performed a CRISPR/Cas9 screen targeting 6237 single sites in the human genome. Next, we developed a data-driven framework to design CRISPR/Cas9 sgRNA (crRNA) pools targeting specific chromatin features. The chromatin context was defined using ChromHMM states, Lamin-B1 DAM-iD, DNAseI hypersensitivity, and RNA-sequencing data. We computationally designed 6 distinct crRNA pools, each containing 10 crRNAs targeting the same chromatin state. We show that the toxicity of a DSB is highly similar across the different ChromHMM states. Rather, we find that the major determinants of toxicity of a sgRNA are cutting efficiency and off-target effects. Thus, chromatin features have little to no effect on the toxicity of a single CRISPR/Cas9-induced DSB.

IHC-based Ki67 as response biomarker to tamoxifen in breast cancer window trials enrolling premenopausal women.

Window studies are gaining traction to assess (molecular) changes in short timeframes. Decreased tumor cell positivity for the proliferation marker Ki67 is often used as a proxy for treatment response. Immunohistochemistry (IHC)-based Ki67 on tissue from neo-adjuvant trials was previously reported to be predictive for long-term response to endocrine therapy for breast cancer in postmenopausal women, but none of these trials enrolled premenopausal women. Nonetheless, the marker is being used on this subpopulation. We compared pathologist assessed IHC-based Ki67 in samples from pre- and postmenopausal women in a neo-adjuvant, endocrine therapy focused trial (NCT00738777), randomized between tamoxifen, anastrozole, or fulvestrant. These results were compared with (1) IHC-based Ki67 scoring by AI, (2) mitotic figures, (3) mRNA-based Ki67, (4) five independent gene expression signatures capturing proliferation, and (5) blood levels for tamoxifen and its metabolites as well as estradiol. Upon tamoxifen, IHC-based Ki67 levels were decreased in both pre- and postmenopausal breast cancer patients, which was confirmed using mRNA-based cell proliferation markers. The magnitude of decrease of Ki67 IHC was smaller in pre- versus postmenopausal women. We found a direct relationship between post-treatment estradiol levels and the magnitude of the Ki67 decrease in tumors. These data suggest IHC-based Ki67 may be an appropriate biomarker for tamoxifen response in premenopausal breast cancer patients, but anti-proliferative effect size depends on estradiol levels.

Targetable ERBB2 mutation status is an independent marker of adverse prognosis in estrogen receptor positive, ERBB2 non-amplified primary lobular breast carcinoma: a retrospective in silico analysis of public datasets.

BACKGROUND: Invasive lobular carcinoma (ILC) accounts for 10-15% of primary breast cancers and is typically estrogen receptor alpha positive (ER+) and ERBB2 non-amplified. Somatic mutations in ERBB2/3 are emerging as a tractable mechanism underlying enhanced human epidermal growth factor 2 (HER2) activity. We tested the hypothesis that therapeutically targetable ERBB2/3 mutations in primary ILC of the breast associate with poor survival outcome in large public datasets. METHODS: We performed in silico comparison of ERBB2 non-amplified cases of ER+ stage I-III primary ILC (N = 279) and invasive ductal carcinoma (IDC, N = 1301) using METABRIC, TCGA, and MSK-IMPACT information. Activating mutations amenable to HER2-directed therapy with neratinib were identified using existing functional data from in vitro cell line and xenograft experiments. Multivariate analysis of 10-year overall survival (OS) with tumor size, grade, and lymph node status was performed using a Cox regression model. Differential gene expression analyses by ERBB2 mutation and amplification status was performed using weighted average differences and an in silico model of response to neratinib derived from breast cancer cell lines. RESULTS: ILC tumors comprised 17.7% of all cases in the dataset but accounted for 47.1% of ERBB2-mutated cases. Mutations in ERBB2 were enriched in ILC vs. IDC cases (5.7%, N = 16 vs. 1.4%, N = 18, p < 0.0001) and clustered in the tyrosine kinase domain of HER2. ERBB3 mutations were not enriched in ILC (1.1%, N = 3 vs. 1.8%, N = 23; p = 0.604). Median OS for patients with ERBB2-mutant ILC tumors was 66 months vs. 211 months for ERBB2 wild-type (p = 0.0001), and 159 vs. 166 months (p = 0.733) for IDC tumors. Targetable ERBB2 mutational status was an independent prognostic marker of 10-year OS-but only in ILC (hazard ratio, HR = 3.7, 95% CI 1.2-11.0; p = 0.021). Findings were validated using a novel ERBB2 mutation gene enrichment score (HR for 10-year OS in ILC = 2.3, 95% CI 1.04-5.05; p = 0.040). CONCLUSIONS: Targetable ERBB2 mutations are enriched in primary ILC and their detection represents an actionable strategy with the potential to improve patient outcomes. Biomarker-led clinical trials of adjuvant HER-targeted therapy are warranted for patients with ERBB2-mutated primary ILC.

IL6/STAT3 Signaling Hijacks Estrogen Receptor α Enhancers to Drive Breast Cancer Metastasis.

The cytokine interleukin-6 (IL6) and its downstream effector STAT3 constitute a key oncogenic pathway, which has been thought to be functionally connected to estrogen receptor α (ER) in breast cancer. We demonstrate that IL6/STAT3 signaling drives metastasis in ER+ breast cancer independent of ER. STAT3 hijacks a subset of ER enhancers to drive a distinct transcriptional program. Although these enhancers are shared by both STAT3 and ER, IL6/STAT3 activity is refractory to standard ER-targeted therapies. Instead, inhibition of STAT3 activity using the JAK inhibitor ruxolitinib decreases breast cancer invasion in vivo. Therefore, IL6/STAT3 and ER oncogenic pathways are functionally decoupled, highlighting the potential of IL6/STAT3-targeted therapies in ER+ breast cancer.

The NF-κB Pathway Promotes Tamoxifen Tolerance and Disease Recurrence in Estrogen Receptor-Positive Breast Cancers.

The purpose of this study was to identify critical pathways promoting survival of tamoxifen-tolerant, estrogen receptor α positive (ER+) breast cancer cells, which contribute to therapy resistance and disease recurrence. Gene expression profiling and pathway analysis were performed in ER+ breast tumors of patients before and after neoadjuvant tamoxifen treatment and demonstrated activation of the NF-κB pathway and an enrichment of epithelial-to mesenchymal transition (EMT)/stemness features. Exposure of ER+ breast cancer cell lines to tamoxifen, in vitro and in vivo, gives rise to a tamoxifen-tolerant population with similar NF-κB activity and EMT/stemness characteristics. Small-molecule inhibitors and CRISPR/Cas9 knockout were used to assess the role of the NF-κB pathway and demonstrated that survival of tamoxifen-tolerant cells requires NF-κB activity. Moreover, this pathway was essential for tumor recurrence following tamoxifen withdrawal. These findings establish that elevated NF-κB activity is observed in breast cancer cell lines under selective pressure with tamoxifen in vitro and in vivo, as well as in patient tumors treated with neoadjuvant tamoxifen therapy. This pathway is essential for survival and regrowth of tamoxifen-tolerant cells, and, as such, NF-κB inhibition offers a promising approach to prevent recurrence of ER+ tumors following tamoxifen exposure. IMPLICATIONS: Understanding initial changes that enable survival of tamoxifen-tolerant cells, as mediated by NF-κB pathway, may translate into therapeutic interventions to prevent resistance and relapse, which remain major causes of breast cancer lethality.

Age-correlated protein and transcript expression in breast cancer and normal breast tissues is dominated by host endocrine effects.

The magnitude and scope of intrinsic age-correlated and host endocrine age-correlated gene expression in breast cancer is not well understood. From age-correlated gene expression in 3,071 breast cancer transcriptomes and epithelial protein expression of 42 markers in 5,001 breast cancers and 537 normal breast tissues, we identified a majority of age-correlated genes as putatively regulated by age-dependent estrogen signaling. Surprisingly, these included genes encoding the chromatin modifier EZH2 (which had a negative age correlation) and associated H3K27me3 (which had an inverse, positive age correlation). Among The Cancer Genome Atlas lung, thyroid, kidney and prostate transcriptomes, the largest overlap with breast cancer in age-correlated transcripts was lung cancer, for which about one-third of overlapping age-correlated transcripts appeared to be estrogen regulated. Age-quartile-stratified outcomes analysis of 3,500 breast cancers using EZH2, H3K27me3, FOXA1 and BCL2 proteins revealed distinct age-related prognostic significance. Age correlation in gene expression may thus be an important factor in ER, EZH2, H3K27me3 and other biomarker assessment and treatment strategies.

Characterizing steroid hormone receptor chromatin binding landscapes in male and female breast cancer.

Male breast cancer (MBC) is rare and poorly characterized. Like the female counterpart, most MBCs are hormonally driven, but relapse after hormonal treatment is also noted. The pan-hormonal action of steroid hormonal receptors, including estrogen receptor alpha (ERα), androgen receptor (AR), progesterone receptor (PR), and glucocorticoid receptor (GR) in this understudied tumor type remains wholly unexamined. This study reveals genomic cross-talk of steroid hormone receptor action and interplay in human tumors, here in the context of MBC, in relation to the female disease and patient outcome. Here we report the characterization of human breast tumors of both genders for cistromic make-up of hormonal regulation in human tumors, revealing genome-wide chromatin binding landscapes of ERα, AR, PR, GR, FOXA1, and GATA3 and enhancer-enriched histone mark H3K4me1. We integrate these data with transcriptomics to reveal gender-selective and genomic location-specific hormone receptor actions, which associate with survival in MBC patients.